EISSN 2319 — 7765 


ISSN 2319 - 7757 


RESEARCH 


Indian Journal of 


Engineerin 


Application of design of experiments (DOE) to simulate selective 
laser melting process for optimum temperature distribution 


Vijay S'*, Gopala Krishna A?, Venkateswara Rao M' 
1. Dept. of Mech. Engg. Bapatla Engg College, Bapatla-522101, India 
2. Associate Professor, JNTU College of Engg. Kakinada, India 


*Corresponding Author: Dept. of Mech. Engg. Bapatla Engg College, Bapatla-522101, India, E-mail: vizzu515@yahoo.co.in, Ph: 
+91 9492466076 


Received 16 September; accepted 13 September; published online 01 November; printed 16 November 2012 


Indian Journal of Engineering, Volume 1, Number 1, November 2012 


ABSTRACT 


The emerging technology in manufacturing is the Rapid Prototyping, by which one can fabricate the functional testing parts within less time. Rapid 
Prototyping (RP) includes different techniques viz., Stereolithography (SLA), Selective Laser Sintering (SLS), Fused Deposition Modeling (FDM), 3D 
Printing etc.. By the advancement of laser technology a new method called selective laser melting is evolved to fabricate metal parts directly from the 
metal powder. The difference between selective laser sintering(SLS) and selective laser melting(SLM) is that with SLS the powders are combined by a 
binding admixture, while with SLM the laser fully melts the powders, producing near full density parts in one step. Temperature distribution while 
fabricating parts plays a vital role both in physical properties and mechanical properties. Model simulation using Analysis packages is a new trend in 
research to minimize time and experimental cost. In this paper a powerful statistical tool called Design of Experiments is used to simulate the Selective 
Laser Melting Process to find the dependency of various parameters like Laser Power, Scan Velocity, scan interval for optimum temperature 
distribution.9OW-7Ni-3Fe powder thermo-physical properties are used for simulation using ANSYS. Three levels for each process parameters were 
selected based on the literature and L9 orthogonal array design is used to investigate the most significant parameters of the process for maximum 
temperature distribution. 

Key words: Selective Laser Melting- Simulation model- Design of Experiments- Temperature Distribution. 


1. INTRODUCTION 


Design of Experiments (DOE) is one of the most powerful statistical tool to investigate the optimum solutions with less effort both in time and cost. 
Application of this method is vast is agricultural, production quality research. Comparatively less number of applications is observed in engineering 
research. R.Anitha et al. [1] used taguchi technique to investigate the critical parameters influencing the quality of prototypes of Fused Deposition 
Modelling. H.J.Yang el al. [2] finds the shrinkage compensation of the SLS process by using the Taguchi method. K.Chockalingam et al. [3] used DOE 
technique to optimize the strength of parts made by Stereolithography process.Mircea Ancau and Cristian Caizar [13] used Pareto-optimal set for 
multicriterial optimization of rapid prototyping processes. Basic steps in Design of Experiments are as follows [12] 

(a) Recognition of and statement of the problem 

(b) Choice of factors, levels and ranges 

(c) Selection of response variable 

(d) Choice of experimental design 

(e) Performing the experiments 

(f) Statistical analysis of the data 

(g)Conclusions and suggestions 

Model simulation is a well known technique to study the behavior of physical systems without disturbing the same. D.E. Karalekas et al. [4] used 

photo elastic models of Stereolithography processed parts for stress analysis investigations. R.K. Chin et al.[5] proposed a thermo mechanical model for 
molten metal droplet solidification. Ruidi Li et al. [7] proposed a 3D transient thermal finite model to p 
simulate the temperature distribution of Selective laser melting process. D.Q. Zhang et al. [6] used 3D 
finite element model to simulate select laser melting of W-Ni-Fe powders. 
Simulation is imitation to the real world problem. Once the simulation model is developed which behaves 
as the physical model, that model can be used to study the response of the physical model for various 
controllable factors. Where as in Design of Experiments, most significant factor on the process can be 
found by doing the experiments on the physical system itself. In this paper a new approach is proposed to 
combine the computer simulation model and Design of Experiment to optimize the temperature distribution 
of Selective Laser Melting (SLM) so that merits of both the techniques are acquired in investigation. 


2. FORMULATION OF HEAT TRANSFER MODEL FOR SLM 


The model of SLM processing based on FE method is shown in Fig.1. The substrate is fixed on the inner 
of forming chamber by four symmetrical screws and is restricted to moving only in the direction of z- axis. 
The model consists of two different parts: (a) a metal substrate (underneath the model); (b) a powder bed 
(upside the model), with the dimensions of 2 x 3 x 1.5 mm and 1 x 2 x 0.05 mm. In order to obtain exact 
calculation as well as sufficient calculation efficiency, the substrate is meshed by solid 90 as a free grid 
while the powder bed is meshed by solid 70 elements with the dimension of 0.05 X 0.05 X 0.05 mm. Figure 1 

ANSYS 12.0 software is used for transient heat transfer analysis of the model described. The 3D heat SLM Model used for analysis 
transfer equation is [9] 
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Where C,, = [N][N]' dv is the heat capacity matrix, [N] is the shape function matrix, [K ] = pc[B] [B]' dV is the heat 
F pc T 
V V 


conduction matrix, {T} and T are, respectively, nodal temperature vector and nodal temperature rate vector, and {Q} is heat flux vector. 


The assumptions used in this analysis are 
° The heat convection transfer only happens between outer boundary and atmosphere 
° The chemical reaction, fluid, and characteristic in the tiny melting pool was neglected 
° The laser power density submits the distribution of Gaussian heat source. 


Table1 Thermal properties of W-Ni-Fe metal powder 


Temperature (T CK) 200 600 800 1000 2000 3500 
Heat conductivity (A Wm! K") 141.9 112.18 110 103 95 83 
Specific heat (Cp Jkg! K") 137.41 144.41 147.84 151.74 | 183.53 | 294.36 


3. MATERIAL PROPERTIES USED FOR THE ANALYSIS 


Investigations of K. Dai and L. Shaw [8] show that there is close relationship between the temperature dependent thermal physical parameters and 
temperature field distribution, especially heat conduction and effective thermal conductivity. There is good mutual solubility between Ni and Fe which lead 
to a complete solid solution. In addition, the large solubility of W in Ni and Fe made it easier to form y- substrate phase with a relative lower melting 
temperature. A homogeneous alloy characterization can be found with a special W-Ni-Fe mixture powder according to W: Ni: Fe equal to 90:7:3 by 
weight. Using the theory of calculating parameters of mixture powder, the thermal physical parameters of composite metal powder are calculated using 
Eq.2 [9] and are as shown in Table 1. 


2 
Cp, =3R| 1- % ~ |+ 2bT + 4dT* — (2) 
20T 
Where a. is Debye characteristic temperature, when using Cp [J/Kg. K] and T(K) to express, Eq.2 can be rewritten as below[9]: 
C, =0.1357676x10*| 1- = + (9.1163 x10 )T + (2.313154 x10°)T? @) 


Investigations of D.Q. Zhang and Q.Z Cai, et al., [9] shows the coefficient of heat conductivity takes great effect in the forming mechanism. It can be 
calculated using Eq. 4[9], if there is no deformation between the metal powders with the shapes of balls: 


k. = Px, 2 1 k, k, 
== Jl-¢@ i +J1-@x 7 i In ; as — (4) 


g g 1-2) T=: 8g 8g 
s S 


Where kg and ks stands for heat conductivity of the environment gas and solid material respectively. © is the material porosity about 0.477. kr is the heat 
conductivity coefficient introduced by radiation of the metal powder 


3.1. Moving Gaussian heat source model 
In SLM technique, the laser energy can be regarded as the form of heat flow density imported into the powder bed which obeys the Gaussian heat source 
distribution [9] 


2AP 2r? 


q= CED) |, ee (5) 
0 oo 


Where w is equivalent radius of laser beam, the distance is from the center of the laser beam area to the point with heat flow density mitigated to the 1/e?. 
A is powder bed absorptivity to laser beam; P is the laser power, the distance from powder bed to the center of the laser beam area can be written as [9] 


Pax e-Bay (6) 


3.2. Treatment of the boundary conditions and initial conditions 
For the domain D, the initial condition is 


T(x, y,z,0)=T, for(x,y,z)¢D 
The natural boundary condition can be defined by 


a, q+ WO -1,)+08(T" =7")=0 
nN 
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On the boundary, S represents those surfaces that are subject to radiation and convection and impose heat fluxes. Where Ae is thermal conductivity 
normal to S; h is heat transfer coefficient for convection; o is Stefan-Boltzmann constant for radiation; ¢ is emissivity and To is the ambient temperature. 


4. DESIGN OF EXPERIMENTS TO SLM 


The different controllable factors of SLM process are velocity of the laser beam v (mm/s), power 


f Table 4 Repose for L9 orthogonal arra 
Table 2 levels of different factors in SLM used for laser beam P (W), Scanning 
interval S (mm), 
Preheat temperature Tp, Ambient 


A. Power(W) 70 85 100 temperature To etc., According to the 4 oe aan 
B. Velocity(mm/s) 20 80 140 investigations of D.Q.Zhang et al., [9] 3 492 4500 
C. Scanning Interval (mm) | 0.05 | 0.075 | 0.1 and LONG Ri-sheng et al., [11] the 
D. Pre heat Temp. (°C) 100 | 150 200 temperature distribution is mostly 4 S77 9000 
influenced by Power, velocity, 5 818 10000 
Table 3 L9 orthogonal arra scanning interval and ambient 6 567 4100 
“sud | run | Feett | Fanz | Fees | raat temperature. Three levels for each 7 912 8000 
factor are selected as in Table 2, in 8 877 8000 
2 1 a a c3 eh which the ranges are selected 9 393 3600 
! 2 al a C1 De according to the literature. 
; 3 fe = & = The objective function of the model is to maximize the temperature distribution of the 
6 5 AD B3 C1 D2 temperature in SLM process. This can be described by the temperature gradient of the moving 
9 6 a3 B3 C2 Di heat source in the plane from the point of application of the heat. If the temperature gradient is 
8 7 a3 Bo C1 D3 more, means that isothermal curves will move farther from the source point. With four factors 
7 8 a3 Bi C3 D2 and three levels of full factorial design includes 64 runs which is time taking and expensive 
2 9 re BD Co D2 in practical. Taguchi’s orthogonal array of L9 series provides more accurate assessment for 


finding the most effective main factors on the response with less number of runs. The 
experimental design for L9 orthogonal array as shown in Table 3 is used for simulation. In the above Table 2 indicates low level, 2 for medium and 3 for 
high levels of the factors. Random numbers for the simulation run is generated through the statistical technique and is implemented for the runs. In the 
first run 85 W of power, velocity of laser beam as 80 m/s, scan interval of 0.1 mm and pre heat temperature of 100°C is used as the input variables with 
material properties and element type unchanged in the simulation using ANSY transient heat transfer. Simulation is carried out for two turns of deposition. 
The response for each of the simulation is temperature distribution in terms of temperature difference and temperature gradient and is shown in Table 4. 


5. RESULTS AND DISCUSSION 


APDL programme is developed to investigate the temperature distribution and gradient for all the runs. Fig. 2 shows the temperature distribution of run 6 
in which power of laser is 100 W, scan velocity of 140 mm/s, scan interval of 0.075 mm and pre heat temperature of 100°C, at the time of laser beam 
reaches the end of one 
complete turn. Fig 3(a), 
(b), (c) shows the 
temperature variation of 
first row middle node 
while the laser beam 


TEMP. DISTRIBUTION FOR MIDDLE NODE 


3500 4 moves for one complete 

3000 | turn. Here the three 
Se —TEST2 graphs are drawn based 
© 2500 + 2 
a —=test3|| on the scan velocity. Fig 
fi 2000 —rstell 3a), (b) and (c) are for 


1500 4 beam velocities of 20 
mm/s, 80 mm/s and 140 
mm/s respectively. 

For knowing the 
most influencing factors 


0 


0.00E+00 5.00601 1.006+00 4.50800 2.006100 2.5000 for the optimum 

TIME (sec) temperature distribution, 

statistical technique 

sa iets Figure 3(a) called Analysis _ of 
gure ——_—— Temp. variation of middle node for run 2,3 and 8 Variance (ANOVA) Is 
Temperature distribution of run 6 used which is a tool in 


statistical software 
Design Expert (DX-8.0). 
Result analysis is given 
3500 3500 in the Table 5, from this 
power of the laser beam 
(A) and scan interval (C) 
ean ay are identified as most 


TEMP DISTRIBUTION FOR MIDDLE NODE TEMP. DISTRIBUTION FOR MIDDLE NODE 


2000 —rTst1 8 2000 —TesT4]/ significant factors for the 
= 1500 — © 500 —TSTé!|_ optimum response. The 
. e significant factors can 
me 1000 also be identified using 
am 500 half normal graph as 

0 0 shown in Fig 4. In the 
0.00E+00 1.00E-02 2.00E-02 3.00E-02 4.00E-02 5.00E-02 6.00E-02 0.00E+0 5.00E 1.00E 1.50E 2.00E 2.50E 3.00E 3.50E above graph factors A 
TIME(sec) 0 0302 PimE(se8F = FRB and C are away from 

the normal line which 

Figure 3(b) Figure 3(c) means that a small 
Temp. variation of middle node for run 1,7 and 9 Temp. variation of middle node for run 4,5 and 6 Naligiion:. ot ites 


factors will influence the 
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response more. The outcome of the Design 
of experiments is the response equation 
which is in terms of factors with its 
corresponding coefficients. From this 
response equation one can estimate the 
predicted values for any level of factors 
within the range. In finalizing the most 
significant Design Expert software is 
performing other tests like residual test, Box- 
Cox plot for power transformations, variation 
of one factor at a time etc., Fig. 5 shows the 
normality of the residual of the model, which 
shows that the residual is propagating on 
either sides of the normal value of the model. 
Fig. 6 shows the graph of residual versus 
run, the residuals are scattered about the 


Figure 4 Figure 5 exact value of the model, which is a 
ar ; Normal plot of residuals for all the runs necessary condition for the runs in DOE. Fig. 
Half normal plot for significant main effects p 7 show the graph of predicted versus actual 


response of the model. The model selected 
for analyzing the data using Design Expert is 
logio for accommodating the large of range 
of data. The graph shows that the predicted 
values are very close to the actual values. 


6. CONCLUSION 


A powerful statistical tool called Design of 
Experiments is used to investigate the 
significant factors for Sintered Laser Melting 
process. Simulation of the process is carried 
out using FEA software called ANSYS and 
the results were analyzed using Design 
Expert. From the analysis it is observed that 
for optimum temperature distribution, power 
Figure 6 Figure 7 of the laser (A) and scanning interval (B) are 
Checking the scatterings of residuals Predicted vs actual values of the response model most significant factors. Care must be taken 
while changing these variables because a 


small variation of in these factors level will affect the response more. 
Table 5 ANOVA for L9 array 
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